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Abstract 

We present an extension of the so-called cumulant crossing method which is 

used for determination of critical point in Monte Carlo simulations. The new 

method uses linear combination of several different order-parameter moments 

and almost eliminates the systematic deviation of crossing points from the 

true critical point. The performance of the method is tested by applying it 

to the 2D Ising model. 

PACS numbers: 68.35.Rh, 52.65.Pp 
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Typeset using REVTgX 



The phenomenological renormalization method [|TJ has been widely used in Monte-Carlo 
simulations of critical phenomena. It make use of the property of reduced order-parameter 
moments such as 

UP =< M n > L /< M 2 > n L /2 (1) 

where M is the order-parameter and < • ■ ■ >l denotes the thermal average in finite system 
with linear size L. In the vicinity of the critical point, they obey following scaling form [|TJ: 

Ut ] ~ f n ,o( tLl/u ) + L-^f n>1 (tL^) + L'^f rii2 (tL 1 ^) + ■■■ (2) 

where t is the temperature deviation from the critical point and v is the critical exponent 
related to the correlation length. The correction-to-scaling and correction-to-correction-to- 
scaling exponent etc. |2| are denoted by u>i(0 < lo\ < u 2 < ■ ■ •) and f n ,i denotes associated 
scaling functions. Two quantities, and uf%; , intersect at the following temperature: 

t ~ dL- 1 ^ 1 + C 2 L' 1/v - UJ2 + ■■■ (3) 

r (l-b-^)fnAO) U] 

1 ~ (bv» - i)/; (o) l4j 

One can extrapolate the critical point by plotting the crossing temperature versus L~ l / U - U ~t 
. However, this method requires a knowledge about the exponents v and ui\ before the 
estimation of the critical point. 

We present an extended method which reduces the systematic deviation of the crossing 
points from the true critical point and does not require any finite-size scaling. The key 
idea is that one can eliminate the correction-to-scaling terms by using linear combination of 
several different moments. 

V L = a&f + a e U^ + ■■■ (5) 
It obeys the following scaling form: 

V L ~ F Q {tL 1/u ) + L-" 1 F 1 (tL 1 ' u ) + L^ 2 F 2 {tL 1/u ) + ■■■ (6) 



where Fi(x) = 04/4^ (x) + a$f§ t i{x) + • • • is the scaling function of Vl- One can adjust the 
coefficients a n so that equations i*i(0) = 0, F 2 (0) = 0, ■ ■ • holds, which leads to C\ = 0, C 2 = 
0, ■ • • in the equation^). For example, consider the following simple case: 

Vl = a,uf + a 6 C/f (7) 

The systematic deviation can be reduced by adjusting the coefficients a n ; the deviation is 
approximately proportional to: 

04/4,1(0) + Q 6/6,i(0) I -i/^ 1 (8) 
a4/l,o(0)+a 6 /4o(0) 1 ' 

Moreover, one can control the direction of the deviation by adjusting the coefficients a n and 
determine the upper and lower bound of the critical point. 

However, the estimation of appropriate coefficients a n is a non-trivial part of this method. 
We show an example of the way of estimation in the following. 

If we plot uf^ versus using the values at the critical point, all points are approxi- 
mately on a straight line. 

(L#>(0), C/f(0)) ~ (/ 4 , (0), / 6i0 (0)) + L-^(/4,i(0), / 6 ,i(0)) (9) 

The slope of the line gives the ratio of two constants /4.1(C)) and /6,i(0). The slope at 
the critical point can be approximated by the value at a temperature where the following 
condition is satisfied: 



(10) 



where b 7^ c 7^ 1 is some constant. 

We have applied this method to Monte Carlo data of the 2D Ising model. Simulations 
were done at the critical coupling (K c = log(l + v2)) using Swendsen-Wang algorithm 
0] and 10 5 observations were made after thermalization. Positions of crossing points were 
calculated with reweighting techniques ||. Statistical errors were calculated with jackknife 
procedure. 



Order parameter moments were calculated using cluster size distribution via the following 
equations: 

< M 6 >l = 16 < S 6 > L -30 < S 4 > L < S 2 > L +15 < S 2 > 3 L (11) 

< M A >l = -2 < S 4 > L +3 < S 2 > 2 L (12) 

< M 2 >l =<S 2 > L (13) 

S n = L- 2n ^m] (14) 
j 

where rrij is the size of j'th cluster in the cluster generation step of Swendsen-Wang algo- 
rithm. 

Fig.l shows the size-dependence of the deviations of crossing points in the two cases of 
using and Vl = 517^ + uf^ with b = 1/2. One can see a significant reduction of 
deviation by the use of Vl- The residual deviations of the case of Vl may be the effect 
of higher order corrections. Further extension for eliminating higher order corrections is 
straightforward, but it needs observation of higher order moments of the order parameter, 
which is more subjected to statistical errors. 

This work was supported by a Grant-in- Aid for Scientific Research by the Ministry of 
Education, Science and Culture. 
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Figure Caption 

Fig.l. Plot of crossing points in the case of using conventional (U^) and new {V£) 
parameter. 
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